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Abstract 

We consider a kinetic model of self-propelled particles with alignment interaction 
and with precession about the alignment direction. We derive a hydrodynamic 
system for the local density and velocity orientation of the particles. The system 
consists of the conservative equation for the local density and a non-conservative 
equation for the orientation. First, we assume that the alignment interaction is 
purely local and derive a first order system. However, we show that this system 
may lose its hyperbolicity. Under the assumption of weakly non-local interaction, we 
derive diffusive corrections to the first order system which lead to the combination 
of a heat flow of the harmonic map and Landau-Lifschitz-Gilbert dynamics. In 
the particular case of zero self-propelling speed, the resulting model reduces to 
the phenomenological Landau-Lifschitz-Gilbert equations. Therefore the present 
theory provides a kinetic formulation of classical micromagnetization models and 
spin dynamics. 
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1 Introduction 

The setting of this paper is the same as in [10] and considers a kinetic model of self- 
propelled particles in three dimensions and its hydrodynamic limit. The particles are 
supposed to move with the same constant speed and their velocity orientation (which are 
vectors of the two-dimensional sphere S 2 ) tend to align to the local average orientation, 
with some noise. This model has been proposed as a kinetic formulation of the Vicsek 
particle model [28]. In this paper, we investigate how the addition of a precession force 
around the local average orientation modifies the resulting kinetic and hydrodynamic 
models. 

The derivation of the hydrodynamic model proceeds like in [10] and relies on both the 
determination of local equilibria in the form of Von-Mises-Fischer (VMF) distributions 
and on the new concept of 'Generalized Collision Invariants' (GCI) which allows for a 
non-conservative closure. In the present paper, we show that the VMF equilibria are 
still equilibrium states of the system with precession and we extend the derivation of the 
GCI's to this system. As in [JO], the resulting hydrodynamics (later on referred to as 
hydrodynamic model for self-alignment with precession) consists of a conservation law for 
the density and a non- conservative equation for the velocity orientation. By contrast to 
[10], the velocity orientation equation contains an additional term which accounts for the 
influence of the precession force and which results in the loss of hyperbolicity of the model 
as soon as the precession force is non-zero. This loss of hyperbolicity occurs when waves 
propagate in the same direction as that of the average velocity direction. 

In order to stabilize the model, we introduce a weak non-locality in the alignment 
interaction: the local average which determines the alignment direction is now taken over 
a certain spatial extension around the particle. If the spatial extension of this local aver- 
age is scaled like the square root of the alignment mean free path e, then, the alignment 
direction possesses a small correction of order e to the local average direction. This cor- 
rection involves the Laplacian of the velocity orientation and results in diffusion terms 
in the macroscopic equation. Due to the precession force, this diffusion term has two 
contribution: the first one corresponds to the heat flow for harmonic maps ; the second 
one is in the form of a Landau-Lifschitz-Gilbert term. In the particular case of zero speed, 
the resulting model reduces to the phenomenological Landau-Lifschitz-Gilbert equations. 
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Therefore the present theory provides a kinetic formulation of classical micromagnetiza- 
tion models and spin dynamics [3]. 

In the non-local alignment interaction, we also include a small (of order e) repulsive 
force which, in the macroscopic model, gives riseto both an additional pressure force and 
a pressure precession term. 

The Vicsek model [28] has been proposed as a model for the interaction of individuals 
among animal societies such as fish schools, bird flocks, herds of mammalians, etc (see also 
[U El El E] ) • It is a particle model (or 'Individual-Based Model' or 'Agent-Based model') 
which consists in an Ordinary Differential Equation system for the particle positions and 
velocities. As the number of particles increases, its computational cost increases more 
than linearly and becomes prohibitive for large particle numbers. The Vicsek model is 
set up in a time-discrete framework. A time-continuous version of the Vicsek model has 
been proposed in [TU] and a kinetic model has been deduced from it on formal grounds. 
The rigorous derivation of this kinetic model has been performed in [3J. In the present 
paper, we will directly start from the kinetic level. 

For the modeling of systems consisting of a large number of individuals, it is more 
efficient to use hydrodynamic models, which describe the system by macroscopic averages. 
Many such models exist [3 [111 EQl EU EH E7] but few of them are rigorously derived 
from particle or kinetic dynamics. The first rigorous derivation of the hydrodynamic 
limit of the Vicsek model is done in [10] (earlier phenomenological derivations can be 
found in [181 EH E5] but do not lead to the same model). The resulting hydrodynamic 
model has also been found in the hydrodynamics of the so-called Persistent Turning 
Walker model [IT] of fish behavior, where particles interact through curvature control. 
Diffusive corrections to the hydrodynamic model have also been computed in [T2] through 
a Chapman- Enskog expansion up to the second order. Other macroscopic models of 
swarming particle systems derived from kinetic theory can be found e.g. in [6] . In 
particular, variants of the kinetic model of [10] and its hydrodynamics have been proposed. 
In [15], the influence of a vision angle and dependency of the alignement frequency upon 
the local density has been investigated. In [161 E] , the modification of the normalization 
constant of the alignment force results in phase transitions from disordered to ordered 
equilibria as the density increases and reaches a threshold. In the spatially homogeneous 
case, these phase transition models are analogs of the Doi-Onsager [T3l [22] and Maier- 
Saupe [19] models for phase transition in polymers. This makes a strong difference in the 
resulting macroscopic models as compared with the ones we are presenting here. 

The organization of the paper is as follows. In section El we set up the kinetic model 
and discuss the equilibria and GCI's. In section [3j the hydrodynamic limit is established. 
The hyperbolicity of the resulting model is investigated in section HI Then, the weakly 
non-local model is derived in section [5j Finally, a conclusion is drawn at section [61 
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2 Setting of the problem and preliminaries 
2.1 The kinetic model 

The starting point is the following kinetic model for a distribution function f(x,v,t), 
where x E R 3 , v E § 2 , t > 0: 

f s + cv . VJ e = ls/ v -[-(P v ±n f e)f + dV v f £ + a(Q f sXv)f £ ] (2.1) 
:= -Q{t) 

:= iCQoCT) + *(/*)), 

where 

Qo(f) = v v -[-(p v ±n f )f + dv v f], 

R(f) = aV v -m f xv)f}. 

Qo(f) is the operator studied in [10]. R(f) is the precession operator about an axis which 
is that of the average particle velocity, c > is the speed of the self-propelled particles. 
We let P v ± = Id — v ® v be the projection operator onto the plane normal to v and 

Pf = / f(v) dv, flf - 



J s2 f(v)vdv\' 



the density and mean velocity direction associated to /. We define the equilibria (Von- 
Mises-Fischer distribution): 

exp((3v ■ Q) 



f s2 exp(/9f • O) dv' 
with P — 1/d. We note the compatibility relation: 

PpF n — Pi tipF n = 

We define 

Mf = Pf F Uf , 
the equilibrium distribution associated to /. 

2.2 Equilibria 

In this section, we investigate the solutions of the equation Q(f) = 0, which are the 
so-called local equilibria of the kinetic model. 
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Proposition 2.1 We have: 

Q(f) = o <^ 3p > 0, 3Q G § 2 , swc/i that f = pF n . 



Proof: We note the following properties: 



/ 

Qo(f)jfdv = -d 



_f_ 

Mi 



[ R(f)-Ldv = -a [ (Q f xv) 



M f dv < 0, 



f(v)dv = 0. 



The first line is proved in [10]. To prove the second line, we remark that V v Mf 
(3P v ±cuMf, (Qf x v) ■ (P v ±flf) = and that the first term can be written as 



a 

— / V, 



(n ' x "»4 



dv = 0. 



Then, if Q(f) = 0, we have 



M f dv. 



This implies that is a constant, and that therefore, / is of the form pF n for a given 
p > 0, Q G § 2 . Conversely, by the compatibility relation, it is obvious that / = pF n 
satisfies Q(f) = 0. ■ 



Remark 2.1 We comment on the connection between this model and rod like polymers. 
In the case of rod-like polymers, the interaction operator is written 

V, ■ [(-V„0 + w x v)f + dV v f], 

where <fi(v) is a potential and w is the precession direction. In general, finding equilibria 
for polymer models is difficult 12^ . Here, we deal with the special case where <fi(v) = —v-w 
which makes this computation possible. In general, the computation is possible if the level 
curves of <fi are invariant under precession rotation. 



2.3 Generalized Collisions Invariants 

We recall the definition of a Generalized Collision invariant (GCI) (see [10]). First, we 
define the collision operator 

Q(f, Q) = V V - [-(P v ±n)f + dV v f + a(Q x v)f] , 

where now Q G § 2 is not necessarily equal to Qf. For given Q, Q(-, Q) is a linear operator 
and we will also consider its formal adjoint Q*(-, Q). Then, we have the 
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Definition 2.2 For given Q £ § 2 , a function ipn(v) is said to be a GCI associated to Q 
iff the following property holds: 



Q) ipndv = 0, V/ such that Qf = ±fl 



We note that the condition Qf = ±fi is a linear condition on / which can equivalently 
be written A ■ J g2 f(v) v dv = 0, for all vectors A such that A ■ Q = 0. On the other hand, 
we have 

.Q)^ndv= [ Q*^ u ,n) fdv. 



'{ijjn,ty = A-v, 



j2 J §2 

Therefore, that ipQ is a GCI is equivalent to 

3A £ M 3 such that A ■ Q = and 



or 



(P v ±£l) ■ V„^ n + dA„?/> c - a(tt x u) • V„Vn = A-v. 

In spherical coordinates (with polar axis Q, latitude 9 £ [0,7r] and longitude y? £ [0, 27r]), 
this equation is written (dropping the index Q): 



sin 9 doip + d 



sin # 



sm l 9 



a sin 9 d^ip 



A\ sin 6 1 cos tp + sin 9 sin 99, (2.2) 



where A = (t4 1 ,t4 2 ,0) are the coordinates of A, in a ortho-normal coordinate basis 
(e 1 ,e 2 ,ft). 

Proposition 2.3 The set Cq of the GCI's associated to Q is a vector space of dimension 
3 spanned by {1, ^ > ^ }> where ip^ k ' is associated to A = e^, for k = 1, 2. Furthermore, 



^ k) (9,p) = 4j[ K) {9)cos<p + ij ( 2 K '{9)smip, k = 1,2 



(*), 



and 



(!) ,/,(2) _ ,/,(!) 



Finally, (ipi,^) satisfies the following elliptic system: 



sin SeVi + d 



sin 9 deip^ + d 



—d e (sm 9 deifj\ 



sin 



sin 2 # 



n 



(i) 



asinfl^ = sin0, (2.3) 



^—rde(sm6doil>. 
sin 6 1 



sin 2 # 



(1) 



+ asin^j 1 ^ = 0. 



(2.4) 



(2) (2) 

(ifti ,ip2 ) satisfies the same system with right hand sides and sin 9 for the first and 
second equation respectively. The elliptic system has a unique zero-average solution, which 
defines ip^ uniquely. 
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Proof: The reduction of (12. 2p to fl 2 . 3 [) . ( 12 .4[) follows from using Fourier transform in the 
ip variable. To solve this system, we construct the complex- valued function ip( k \9) = 



^\0) + iip 2 k \9). ifj( 1 ^ satisfies the complex- valued elliptic problem: 



,(fc) 



sm9d e ^ {1) +d 



sin 9 



d e (sin9 dgtfj 



(i) 



sin 2 9 



(1) 



+ ia sin 9 ij)^ 1 ' = sin 9 



whereas i/j^ satisfies the same equation with right-hand side isin9. Since the imaginary 
part of the operator is skew-adjoint, it does not modify the energy identity. Therefore, 
the same theory as in [10] applies and shows the existence and uniqueness of the solution 
in the space of zero average functions. The remaining part of the statement follows easily. 



3 Hydrodynamic limit 

The following theorem establishes the limit e — > 0. 

Theorem 3.1 We suppose that f e — > f° as smoothly as needed. Then, we have 

f = pFn, (3.1) 
where p = p(x,t) and Q = Q(x,t) satisfy the following system 

d t p + c Cl V x (pQ) = 0, (3.2) 
p(d t Q + cc 2 cos<5 (fi • V X )Q + cc 2 sinSQ x ((fi ■ V x )fi)) + cdP n ±V x p = 0. (3.3) 

The coefficient c\ is defined by 

ci = / F n {v) (v ■Q)dv = — — — . (3.4) 

The coefficients c 2 and 5 are defined as follows. First, define at and bk for k = 1, 2 by 

ak = \l F Q (cos 9) $\e) sin 2 9 d9, fc = 1, 2, (3.5) 

* J0£[O,tt] 

b k =\ ! F Q (cos 9) ^(9) cos 9 sin 2 9 d9, fc = l,2, (3.6) 

where ip^ are the GCFs found in proposition ^. 3[ Then, introducing (p a -,9 a ), (pb,9 b ) the 
polar coordinates of the vectors (01,02) and (61,62)? 

a x + ia 2 = p a e i6a , 61 + ib 2 = p b e ie », 

c 2 and 5 are defined by 

02 = -, 5 = 9 b - 9 a . (3.7) 

Pa 
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Remark 3.1 Model A3.2fy . A3. 3\) will be referred to as the hydrodynamic model for self- 
alignment interactions with precession. 

Remark 3.2 We note the following identities: 

(ft • V x )ft = (V a x ft) x ft, 

ft X ((Va, X ft) X ft) = P n ±{V x X ft). 

Remark 3.3 The reason for keeping the dependence upon the speed c explicit is that we 
will later investigate the case c = 0. 

Proof. Letting e — > in (12.11) and using proposition 12. 1\ we get that Q{f £ ) — > = Q(f°), 
and that eq. (13.1 j) is satisfied. Now, we look for the equations satisfied by (p, ft) and for 
this purpose, we use the GCI. 

According to Proposition ^. 3[ 1 is a classical collision invariant. Therefore, multiplying 
(12. ip by 1 and letting e — > 0, we get the mass conservation equation (13.21) . 

Now, we successively multiply (12. ip by the GCI's ip^ k \ k = 1, 2 associated to the polar 
vector Qfe. From the fact that i/j^ E ^n /£ , we get that 

/ Q{f e )^ k Uv= [ Q(f,n f .)^ k Uv = 0, k = l,2, 
Js 2 Js 2 

thanks to definition 12.21 We deduce that 

/ T{f e )i(>Wdv = 0, fc = l,2, 
Jn 2 

where Tf = d t f + cv ■ V x f is the transport operator. Letting e — > 0, we deduce that 

T {k) := [ T(pF n )^ {k) dv = 0, k = 1, 2. (3.8) 
Js 2 

We note that 

=^(0) cos<^ + ^ 2 (0) siny?, (3.9) 
^( 2 ) = -^ 2 (^)cos^ + ^i(^)sin^, (3.10) 

where we define ipi = ipjp and ip2 = for simplicity. 
A simple computation [10J shows that 

T(pF n ) = F n {d tP + c(v ■ V x )p + 0p(v ■ (d t + cv ■ V,)ft)} . 

We decompose 

V = V± + V\\, V ± = Pn±V, V\\ = (v ■ ft)ft. 

In spherical coordinate system, we have 

v± = (sin 8 cos ip, sin 6* sin </?, 0) T , = (0,0, cos 0). 
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Therefore, v± is an odd degree trigonometric polynomial of <p, while v» is an even de- 
gree such polynomial. We decompose T(pFn) into even and odd degree trigonometric 
polynomials in tp using the above defined decomposition of v. We get: 

T(pF n )=T e + T , (3.11) 

t = F n {cv ± ■ v xP + p P v ± ■ (dt + c(v ■ n)n ■ V x .)fi} , 

where we have used that v\\ ■ {(v± ■ V X )Q) = because \Q\ = 1. We can write 
T = Fq v j_ • D{v ■ Q) — Fq sin 9 (cos tp D\(cos 9) + sin tp D 2 (cos 9)), 

with 

d{v ■ n) = cv x p + (3p(d t + c{v ■ n)n ■ v x )n, (3.12) 

and D = (Di, D 2 , D 3 ) are the coordinates of D in the cartesian coordinate basis associated 
to the definitions of the angles (8, tp). 

Now, inserting (13.111) . (13.121) into (13. 8p . and noting that the odd degree polynomial in 
ip vanish away upon integration with respect to tp e [0, 27r], we get, 

T (1) = / F n (cos 9) (cos tp D l (cos 9) + sin tp D 2 (cos 9) ) 

(ipi (9) cos tp + ip2 (9) sin <p) sin 2 9 d9 dp, 
T (2) = / F n (cos 0) (cos y? D x (cos 0) + sin (p D 2 (cos 0) ) 

</<9e [0,7r], i/?6 [0,27r] 

(—ifj 2 {9) cosp + ipi(9) sin 99) sin 2 6* ci^ d<p. 

Performing the tp integration leads to: 

T (1) = -/ F n (cos0) (Z>i(cos0)V>i(0) + L> 2 (cos 0)^(0)) sin 2 0d0, (3.13) 

2 j6»G[0,7r] 

T (2) = -/ F n (cos0) (-£>i(cos0)^ 2 (0) + /^(cos 0)^(0)) sin 2 0d0. (3.14) 

2 J6»G[0,7r] 

Now, introducing the matrices 
and 

M = 

with a fc and 6 fc (Jfe = 1,2) defined at 1ET5|) . (l3~6l) and inserting (l3~12l) into (jOSjl . (l3~T4l 
leads to 

= [AjP^cV^ + /^fi) + /3pc[5]P n x((fi ■ V*)n) = 0. (3.15) 





a 2 


0^ 






& 2 




-a 2 


ai 











: 








y 




^ 









The matrix [a] being a similitude matrix, it is invertible and we can introduce the matrix 

M - W- N - ( Z Z ) • 

and, 

(c 2 i c 22 \ 
-c 22 c 2 i . 
/ 

Multiplying the first two lines of the vector equation (13 . 1 5[) by d[a] -1 , we get: 

p(d t Q + c[C 2 }(n ■ V x )n) + cdP n xV x p = 0. (3.16) 

This is the momentum equation. 

The equation can be transformed using the fact that [b] and consequently [02] are 
similitude matrices. Therefore, we can write 

M = c 2 R 5 , (3.17) 

where c 2 and 5 are given by (13. 7p and 

_ f cos 5 sin 5 \ 

I — sin 5 cos 5 J ' 

is the rotation matrix of angle 5. Then, (13.161) can be equivalently written according to 

(13. 3p . which ends the proof. ■ 



4 Hyperbolicity of the hydrodynamic model 

In this section, we investigate the hyperbolicity of the hydrodynamic model for self- 
alignment interactions with precession. For this purpose, we use the spherical coordinates 
associated to a fixed cartesian basis. In this basis, denoting by 9 G [0, 7r] the latitude and 
ip G [0, 2ir] the longitude, we have 

f2 = (sin 9 cos ip>, sin 9 sin (p, cos 9) T , 

and we let Qe and be the derivatives of Q with respect to 9 and ip. We note that 



\Qg\ = 1, \Qtp\ = sin^. 



We will use the formulas 



v x -n = n e - vj + % ■ v x ^, 

(ttp ■ a) 
sin 2 9 

(q ■ v x )n = {{n ■ v x )9) n e + {{n ■ v x )<p) %, 

fix(n- V X )Q = _ sin . V x )tp) tie, 

sm u 

fit = Q e 6 t + fL tp t , 
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where a is an arbitrary vector. 

As in [10], we use the time rescaling t' = cc\t and introduce 

C2 x2 d A 
a = — , A = — , p = \mp. 

Ci Ci 

With this change of variables and unknowns, system (13.21) . (13. 3p is written (dropping the 
primes on t): 

p t + n- v x .p + av, ■ n = o, 

fit + a(coso"(fi ■ V x )tt + sm5Q x ((fi • V x )fi)) + AP J] iV I| o = 0, 

or, 

p t + fi • V x /3 + A(fi e • V x # + fi v ■ V x ip) = 0, (4.1) 
9 t + a (cos 5 (fi ■ V^)^ - sin 5 sin (fi ■ VJ^) + A fi e - V x /3 = 0, (4.2) 

<p t + a (cos 5 (fi ■ V x )<p + ^ (fi ■ V a )0) + A ^,'7/ = 0. (4.3) 

sin sin & 

In passing, we note a mistake in formula (4.70) of pjj] where the last term should be 
divided by sin 2 9. This mistake does not affect the eigenvalues of the system which are 
correct. 

Introducing 

' P 

U 

this system is written 

U t + A(U)U X + B(U)U y + C(U)U Z = 0, 

in cartesian coordinates x = (x,y,z). Assuming translation invariance along the z direc- 
tion, the system is reduced to 

U t + C(U)U z = 0, 

with 

(cos# —A sin 9 
—A sin 9 a cos S cos 9 —a sin 5 sin# cos# 
a sing cose acoso"cos# 

sin o 

The system is hyperbolic, if and only if the eigenvalues of C(U) are real for all values of 
U. The characteristic polynomial of C(U) is given by 

P(X) = A 3 — (1 + 2a cos 5) cos 9 X 2 + (a (a + 2 cos 6) cos 2 9 — X 2 sin 2 9) X + 

+a (A 2 cos S sin 2 9 — a cos 2 9) cos 9. 

We note that in the case sin 5 = 0, we recover the results of PJJ] where all eigenvalues are 
real. In the present sin 5 ^ case, the roots of P(X) cannot be analytically computed. 
However, they are analytically computable in the two cases 9 = and 9 = tt/2. 
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Case 9 = vr/2. Then, P(X) = X(X - X)(X + A) has real roots X = and X = ±A. 
These roots are the same as in the 5 = case [TO] . 

Case 9 = 0. Then, P(X) = (X - 1)((X - a cos<f) 2 + a 2 sin 2 5) has one real root X = 1 
and two complex conjugate roots X = ae ±lS . Therefore, as soon as sin 5 ^ 0, the system 
loses its hyperbolicity near 9 = 0. 

Due to the loss of hyperbolicity of the hydrodynamic model, we look for diffusive 
corrections by introducing a nonlocal evaluation of the alignment direction. This task is 
performed in the next section. 



5 Taking into account non-locality in the interaction 

Now, we suppose that the Qf vector is replaced by a non-local evaluation denoted by Q v f . 



So, we introduce the following kinetic model: 



ft + cv V x f = -V„ 

e 



-(P v ^ r >)f + dV v f + a(mxv)f 



(5.1) 



where 



\jf + rfr 



j]= K[ l - ^ ) v' f(x', v', t) dv' dx' 



,'/ 





.( 


' \x' — x\\ 


/ 




i(x',v')e 


R3 X § 2 V 


, V ) 



f{x',v',t) dv'dx', 



is a force term. The first term (given by jj) expresses the alignment interaction (like in the 
Vicsek dynamics) while the second term (given by rj) expresses the repulsion interaction. 
We denote: 



K(\£\)dt = k , 



x'£R n 
1 



2n 



K(\Z\)\tfdZ = k, 



x'eR n 



We can always assume that k$ = 1. Repulsive interaction here means that we assume 
(f) > 0. Defining Uf by 



Pf u f = I f(v')v dv' 
>v'e§ 2 
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we have the following Taylor expansion of nj: 

ft} = 0, + if-^-if + o(n 2 ), £ f := P n x(fcA(p /W/ ) - 0V, P/ ). 

p f \Uf\ t 

Inserting this expression into the kinetic equation (12. II) . we get 

ft + cv VJ = ^V v ■ [-(P v ±n f )f + dV v f + a(Q f x v)f] 

2 i 2 

+^^V, ■ [-(P v ,£ f )f + a(£ f x v)f] + oCL). 
e Pf\ u f\ e 



We now let — = 1 and we let e — > 0. Dropping terms of order o(l) or smaller, this 
leads to the following problem: 



ft + cv ■ V x f e = ■ [-(P,iO f )f + dV v f + a{Q f e x v)f] 

+ ^—V v ■ [-(P v± £ fs )f + a{l fe x v)f] (5.2) 

Pfe\Ufe\ 

= -Q{n - L(n, 

with 

Lf = • [-(Pv^f)f + <*{£ f x v)f] . 

Pf\U f \ 

The following theorem establishes the formal e — > limit. 
Theorem 5.1 We suppose that f £ — > f° as smoothly as needed. Then, we have 

f° = pFn, (5.3) 
where p = p(x,t) and ft = Q(x,t) satisfy the following system: 

a t p + cc 1 V x { / M) = 0, (5.4) 
p{d t Vt + cc 2 cos 5 (ft ■ V x )ft + cc 2 sin 5 ft x ((ft ■ V x )ft)} + cdP n ±V x p + 

+— {-(2d + c 2 cos5)P n ±(kc 1 A(pQ) - <pV x p) + 

Cl 

+ (c 2 sin5 - a) (ft x (fcciA(pft) - 0V »)} = 0. (5.5) 

and tae coefficients c\, c 2 and 5 are tae same as in Theroem \3.1[ We restrict ourselves to 
the case 2d + c 2 cos 5 > which is the condition for the system to be stable. 

Remark 5.1 A special case is c = 0, p = 1, <f> — 0. This leads to 

d t VL + k{2d + c 2 cos 5) ft x (ft x Aft) + k(c 2 sin 5 - a) (ft x Aft) = 0, 

which is the classical Landau-Lifschitz- Gilbert equation, provided that 2d + c 2 cos5 > 0. 
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Proof. Since Lf is a divergence in velocity space, the mass conservation equation is 
unaffected by this additional term compared to section [31 The additional contribution of 
Lf to the momentum equation is a term of the form: 

£(*>:= [ L(pF n )^ k) dv, k = 1,2, (5.6) 



at the left hand side of (13.151) . We have 

L(pF n ) = --V„ • i-(P.^£)F n + a(£ x v)F n ] 

Cl 

where 

£ := i pFa := P Q ±(k Cl A(p{l) - </>V x p). 
Therefore, using Green's formula, we get: 



Cl 

Now, we use the formulas: 



£ (fc) = — / \-(P v ±£) + ail x v)\ ■ V v ip {k) F Q dv 
= — / [-£ + a(£xv)]-V v tp {h) F n dv = 

C l Js 2 

- f / [-V^ (fc) + a(v x V^ (fe) )l ^ 
-l Vis 2 



V„<? dv = 2 vg dv, 
Js 2 

/ iy v g)hdv = 2 I vghdv— / (V v h)gdv, 
Js 2 Js 2 Js 2 

(v x V v g)hdv = — / (t> x V v h)gdv, 



for any pair of scalar functions g, h on E> 2 . We get: 

= ~( [ i-2vF n + V v F n -a(vxV v F n )]^ h Uv)-£. 

c l VJsa / 

Now, we note that V„-Fh = ^(P^xfiji^h, which leads to 

£(*) = —( f [-2v + (3(P v ±n)-a(3(vxn)}ij {k) F n dv)-£, 
Cl \>/§ 2 / 

(we used that ?; x = i> x (P„xfi)). We use the decomposition v = v±_ + vu. Since £±Sl 
and Pq±P v ±Q = — (v ■ £l)v±, we can write: 

£(*) = 1 ^jf [-(2 + P(v ■ Q))v± -aP{v± x fi)] 7/; (fe) F n rf^ 
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Introducing the coordinates £ = (£i,£ 2 , 0) in the cartesian basis associated to Q, and 
decomposing v±_ accordingly, we get: 

C {k) = —f [-(2 + /3cos0)(cos^i + smtp£ 2 ) 

— a/3(sin (^£1 — cos (p£ 2 )} -Pb sin # dv 
= — I [(-(2 + /3cos0)£i + a/%)cos<^+ 

+ (-(2 + /3 cos fl)£ 2 - sin <p] ^ {k) F n sin 9 dv. 

Finally, using (ET5]) . fEUO]) . we get: 

= — f [{-{2 + /3cos9)£ 1 + aP£ 2 )cosip+(-(2 + /3cos9)£ 2 -a/3£ 1 )smip} 

(ipi(9) cos ip + ip 2 (9) simp) Fq sin9dv, 
£ (2) = —[ [(-(2 + (3 cos 9)^ + a/3£ 2 ) cos ip + (-(2 + (3 cos 9)£ 2 -a(3£ 1 ) simp] 

C l i§2 

(—^2(6') cos y + sin<p) sin0<fo, 
and performing the integration with respect to tp G [0, 27r], we are led to 

= ^- f [(-(2 + /3cos9)£ 1 + a/3£ 2 )tp 1 (d) + (-^ + ^os9)£ 2 -a/3£ 1 )tp 2 (9)} 
2c\ Jo 



F n sin 2 9 d9, 

£ (2) = ^- f [-{-{2 + /3cos9)£ 1 + a/3£ 2 )iP 2 {9) + {-{2 + /3cos9)£ 2 -af3£ 1 ) 

F n sin 2 9 d9. 

In vector notations, this is written: 





Finally, collecting all terms leads to the following equation: 



£(2) I = I {-(2[A] + f3[B])£ + a(3[A](£ x Q)} 

Cl 



[A]P f ,x(cV :c p + PpdtG) + Ppc[B]P a x((n ■ V X )Q) + 

+- {-(2[A] + f3[B})£ + x ft)} = 0. 

Cl 

Now, multiplying the first two lines of this vector equation by /3[a] _1 , we get 

p(d t Q + c[C 2 ](ft ■ V X )Q) + cdP n ±V x p + — {-{2dld + [C 2 })£ + a(£ x ft)} = 0, 

Cl 
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or, expliciting the expression of t. 

p(d t Q + c[C 2 ](n ■ V X )Q) + cdP n ±V x p + - {— (2dld + [C 2 })P n ±{k Cl A(pQ) - <pV x p)+ 

+ a((fcciA(pfl) - <pV x p) x Q)} = 0, 

Using (I3.17p . we can also write this equation in the form (15. 5p . which ends the proof of 
the theorem. ■ 



6 Conclusion 

In this paper, we have derived the hydrodynamic limit of a kinetic model of self-propelled 
particles with alignment interaction and with precession about the alignment direction. 
We have shown that the resulting system consists of a conservative equation for the local 
density and a non-conservative equation for the orientation. We have observed that this 
system may lose its hyperbolicity. Then, we have provided an extension of the model 
including diffusion terms under the assumption of weakly non-local interaction. In the 
particular case of zero self-propelling speed, we have noted that the resulting model is 
nothing but the classical Landau-Lifschitz-Gilbert equation. Therefore the present theory 
provides one of the very few (if not the only) kinetic justification of the phenomenolog- 
ical Landau-Lifschitz-Gilbert equation. Future works in this direction will include the 
sensitivity analysis of the hydrodynamic model in terms of the modeling parameters, the 
accounting for phase transitions, the influence of a vision angle, and the inclusion of the 
particle interactions with a surrounding fluid in the perspective of modeling active par- 
ticles suspensions. Another research direction consists in investigating the mathematical 
structure of these models, prove well-posedness and investigate its qualitative properties. 
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